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Abstract. Azimuthal correlation functions are used to study jet- and di-jet 
properties as a function of centrality in Au+Au collisions at y/s^7=200 GeV. 
Utilizing a novel technique to decompose the correlation function into a (di-)jet 
and an underlying event, the jet-pair distribution is extracted and compared 
■ to similar results for d+Au collisions obtained at the same collision energy. A 

\ striking similarity is observed between the widths and associated yields of the 

(di-)jet distributions for d+Au and peripheral Au+Au collisions. By contrast, 
the distributions for mid-central Au+Au collisions indicate an increase in the 
di-jet yield with centrality, and a very broad away-side jet having a possible 
minimum at A<f> ~ ir. These features point to significant medium induced 
modification to the away-side jet and are compatible with recent predictions of 
jet-induced "conical flow". 
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1. Introduction 

Reactions between Au ions at the Relativistic Heavy Ion Collider (RHIC), indicate 
the creation of a fireball of nuclear matter having energy density well above that 
required for a de-confined phase of quarks and gluons (QGP) [HE]- The decay 
of this matter results in large azimuthal anisotropics in the particle emission pat- 
terns, suggesting early thermalization and the development of substantial pressure 
gradients which drive the dynamical evolution of the system 0J |SJ [7| 
Hydrodynamic evolution of the fireball is further corroborated by the observation 
of strong radial flow [ 9 and first hints of a long-range emitting source from a recent 
Imaging analysis [ 1101 II 1 1 . Strong indications for hydrodynamic evolution of the 
emitting system implies the production of strongly interacting high energy density 
matter in energetic RHIC collisions [0J El- Indeed, this matter has been ob- 
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served to strongly suppress the yield of hadrons with large transverse momenta [ 
1121 1131 1141 ITS] and to suppress the away-side jet in central Au+Au collisions [I16j. 
It is believed that this suppression results from energy loss of hard-scattered par- 
tons traversing the high energy density matter prior to the formation of hadrons [ 

ei ua ua mi eh ei 

An important open question of great current interest is the influence of the 
parton-medium interactions on jet properties. Such an influence is of paramount 
importance if one wants to use jets as a probe of the properties of strongly interacting 
high energy density matter. Several recent works have outlined a possible influence 
of the coupling between jets and a strongly interacting medium [ 1231 1241 1251 1261 
1271 1281 1291 I3UI I31j . A particularly important proposal is the conjecture that the 
energy deposited in the medium could lead to the creation of a shock wave around 
the propagating parton, thereby creating "conical flow" or "bow waves" analogous 
to a sonic boom in a fluid [ 1771 1281 1291 I3(J| . The experimental observation of such 
conical flow could serve to pin down the sound speed in the nuclear matter created 
at RHIC. 

In order to probe the influence of possible parton-medium interactions on jet 
properties, we use azimuthal angular correlation functions to investigate jet topolo- 
gies and yields in d+Au and and Au+Au collisions. Here, the operational strategy 
is that d+Au measurements provide a good baseline for comparison to the Au+Au 
measurements which are expected to show much stronger modifications to jet prop- 
erties. 



2. Data Analysis 

The analysis presented in this paper uses Au+Au and d+Au data ( v / saT/v=200 GeV) 
provided by RHIC in the second and third running periods (2001,2003), respectively. 
The full PHENIX detector setup is described elsewhere [El]. Charged tracks rel- 
evant to this analysis were reconstructed in the central arms of PHENIX, each of 
which covers 90 degrees in azimuth. Tracking was performed via the drift chamber 
and two layers of multi-wire proportional chambers with pad readout (PC1,PC3)[ 
132) . A combinatorial Hough transform in the track bend plane was used for pattern 
recognition[ 32]. Most conversions, albedo and decays were rejected by requiring 
a confirmation hit within a 2 a matching window in the PC3. Collision centrality 
was determined via cuts in the space of BBC versus ZDC analog response [ 134) . 

The correlation function in relative azimuthal angle between particle pairs, 
A<fi = (0i — 4>2), is defined as the ratio of two distributions 

N cor (A0) 

The foreground distribution N cor , measures coincident particle pairs from the same 
event by pairing particles from a high-py "trigger" bin (2.5 < px < 4.0 GeV/c, 
hereafter labeled A) with associated particles from a lower pt selection (1.0 < pr < 
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2.5 GeV/c, hereafter labeled B). The background distribution N m i X , is generated 
in an analogous way by mixing particle pairs from different events within the same 
multiplicity and vertex class. The Azimuthal acceptance and detector efficiency 
effects cancel in the ratio of foreground to background distributions and the corre- 
lation function yields the probability distribution for detecting correlated particle 
pairs per event within the PHENIX pseudorapidity acceptance (\r]\ < 0.35). Fig.^ 



2.5 < p T Wg < 4.0 GeV/c 4.0 < p T lrig < 6.0 GeV/c 
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Fig. 1. d+Au correlation functions C(A(f>) for several pr selections as indicated. 
The centrality selection is 0-80%. The dashed line represents a double gaussian fit 
to the data. 



shows d+Au correlation functions for two different trigger and associated px selec- 
tions as indicated. The centrality selection is 0-80%. The dashed line represents a 
double Gaussian fit to the data. These d+Au correlation functions exhibit a shape 
reminiscent of what one would expect from di-jet fragmentation. That is, a rela- 
tively narrow near-side peak centered at Acf> = and a somewhat wider away-side 
peak centered at Acf> = ir. Fig. ^ shows that both near-side and away-side peaks 
narrow and are more pronounced for higher ^-selections of trigger and associated 
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particles, respectively. Such a pattern is expected if (di-)jet fragmentation is the 
dominant particle production mechanism at high transverse momenta. It is note- 
worthy that the widths (of near- and away-side jets) and the yields obtained from 
d+Au correlation functions are rather similar to those obtained from p+p collisions. 
Consequently, we ascribe all correlation in d+Au collisions to jets. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 
A (|> (rad) A $ (rad) 



Fig. 2. Correlation functions C(A</>), generated for trigger particles in 2.5 GeV/c< 
Pt <4.0 GeV/c and associated particles in 1.0 GeV/c< p^ <2.5 GeV/c. The bands 
show the harmonic contribution within the systematic uncertainty. The dotted lines 
indicate the value of a a (see text). 



The d+Au correlation functions (cf. Fig. Q are to be compared to the corre- 
lation functions for Au+Au data shown in Fig. |3 The filled squares in the figure 
show correlation functions for AB charged hadron pairs for several indicated cen- 
tralities. It is difficult to overlook the striking similarity between the correlation 
function obtained for the most peripheral collisions (cf. Fig. EJ) and that obtained 
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for d+ Au collisions (cf . Fig. ^) . Both correlation functions show the two narrow 
peaks (located at A(f> = and A<j> — it) characteristic of (di-)jet fragmentation. By 
contrast, the correlation functions for more central Au+Au collisions show strong 
indications for a harmonic component and hints for a broad away-side jet ie. the 
correlation functions show a minimum below Acj) = tt/2. Such a shift away from the 
minimum expected for harmonic contributions, can only come about if the away- 
side jet is significantly broader than that observed in d+Au collisions. An important 
finding that should be stressed here is that the observed characteristics of all of the 
Au+Au correlation functions can be fully accounted for via two contributions to the 
correlation function: (i) a (di-)jet and (ii) a harmonic contribution [ 1161 1351 1361 137j . 

3. Decomposition of jet and harmonic contributions 

Careful investigation of possible modifications to the di-jet distributions in Au+Au 
collisions require access to the jet-pair distribution without the blurring effects of 
the harmonic contributions. Consequently, a reliable procedure for decomposing 
the measured correlation functions into their (di-)jet and harmonic (or flow) contri- 
butions are required. Detailed descriptions of such a procedure are given in Refs. [ 
I3"%1 EH|- We give here only an outline of the main points. 

It can be shown [ I38| that the pair correlations from the combination of flow 
and jet sources is given by 

C AB {A$) = a [C AB {A^)] + J(A0), (2) 

where C^ B ( A(f)) is a harmonic function of effective amplitude V2, 

C^ B {Acp) = [1 + 2v 2 cos2{A4>)]; v 2 = (v A x v B ). (3) 

and J(A</>) is the (di-)jet function. No explicit or implicit assumption is required 
for the functional form of J(A</>). Rearrangement of Eq. Ogives 

J(A0) = C AB (A0) - a C AB {A$). (4) 

Thus, one only requires knowledge about a and V2 to evaluate J(A0). To constrain, 
a we assume that the (di-)jet function has zero yield at the minimum (ZYAM) 
A(fi m i n , in the jet function, i.e. aoC^ B (A<j) m i n ) ~ C AB (A<p m i n ). This fixes the 
value of a Q . 

The V2 value reflects the average anisotropy of the particles from both sources, 
and can be obtained from the single particle distributions relative to the reaction 
plane ipR. However, this step requires that the reaction plane is itself determined by 
a procedure essentially free of non-flow effects. This is accomplished in the present 
analysis by demanding a large (pseudo)rapidity gap (A?; ~ 3 — 3.9) between the 
reaction plane and the particles correlated with it EH It is expected that the v 2 
values so obtained are much less affected by jet contributions [I40|. In addition, the 
reaction plane (at each centrality) and its dispersion correction, and v A and v B were 
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obtained from the same data set used for jet function extraction in order to avoid 
potential biases. Correction for reaction plane dispersion followed the procedures 
outlined in Ref. 

3.1. Simulation tests 

Prior to applying the decomposition method to PHENIX data, its reliability was 
thoroughly tested via extensive Monte Carlo simulations [ |2U ■ These investiga- 
tions included simulation tests which took account of the <fi and the r\ acceptance 
of PHENIX. Representative results from these Monte Carlo investigations are sum- 
marized in Fig. |21 Panel (a) depicts the azimuthal distribution of the simulated 
reaction products with respect to the event plane. The smooth line is a harmonic 
fit to the data to extract v<x. This i>2 determines the amplitude of the harmonic 
component (dashed line in panel (b)) to be subtracted from the correlation func- 
tion (filled circles) shown in panel (b). The solid squares in panel (b) shows the 
ZYAM subtracted jet-pair distribution referenced to a . This distribution is to be 
compared to the input jet pair distribution (solid line) obtained via tagging of the 
jet-particles in the simulation. The rather good agreement shown between input 
and output jet-pair distributions in Fig. 03 serves to confirm the reliability of the 
method. The method is easily generalized to the case in which trigger particle de- 
tection is constrained within a cut angle parallel or perpendicular to the reaction 
plane [EHj. Results from simulations in which such constraints have been applied 
are summarized in panels (c) and (d) of FigEI In this case, the harmonic function is 
determined following the techniques outlined in [I43|. Here again, panels (c) and (d) 
clearly indicate that the input jet function is reproduced in detail. It is noteworthy 
that a wide range of tests for a variety of input jet-pair distributions including those 
that might be expected from conical flow, were made with equally good recovery of 
the input jet-pair distributions [I39|. 

3.2. Decomposition of the measured correlation functions 

The solid bands in each panel of Fig. [21 illustrate the application of the ZYAM 
condition to PHENIX data with the measured values of vi (i>2 = (u^ x Ujf )). The 
dashed lines show the a value obtained for each centrality. Following Eq. 01 the 
jet-pair distribution is obtained at each centrality via subtraction of the harmonic 
contribution from the correlation function. It is straightforward to show that the 
integral of this distribution is related to the average fraction of jet-correlated particle 
pairs per event and hence the conditional per trigger yield [0H1 EES- The ratio of 
the sum of J(A<fi) and the sum of C ( A<j>) (over all bins in A</>) gives the fraction of 
jet-correlated particle pairs per event PF, 



PF = 



EiC(A^) 



(5) 
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Fig. 3. Simulation test of the ZYAM decomposition procedure from [125]. (a) 
simulated azimuthal distribution of particles with respect to the event plane, (b) 
Simulated correlation function with harmonic component and extracted jet-pair 
distribution. The solid line represents the input jet-pair distribution (referenced to 
a ) as obtained from tagging the jet particles in the simulation, (c) and (d) same as 
(b) but for trigger particles constrained in-plane and out-of-plane respectively (see 
text). 



Subsequent multiplication of this fraction by the average number of detected particle 
pairs per event (N^ B ), followed by a division by the product of the detected singles 
rates (N^), (Nf), gives the event averaged jet-pair production in excess of the 
combinatoric pair production. A final product with the efficiency corrected singles 
rate (N^ ff ), for bin B, gives the efficiency corrected pairs per trigger or conditional 
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yield CY [ESI EE], 



4. Results 




A (|) (rad) A § (rad) 



Fig. 4. ZYAM-subtracted jet-pair distributions l/N trig dN/d(A(j)). The 
dashed(solid) histograms indicate the distributions resulting from increas- 
ing(decreasing) V2 — X v^f by one interval of the systematic error. 



The ZYAM-subtracted conditional yield distributions, normalized to the num- 
ber of jet-pairs per trigger particle, are shown in Fig. 0] The distribution for the 
60-90% peripheral event sample (cf. Fig. Ef ) shows the typical (di-)jet shape that is 
familiar from p+p and d+Au collisions at RHIC. It consists of two distinct peaks, 
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one narrow near-side peak centered at A0 = and a broader away-side peak at 
Acf) = ir. It is interesting to trace the evolution of these peaks with collision cen- 
trality, as one expects to make increasing amounts of hot and dense matter in the 
more central collisions. Scanning the mid-central and central jet-pair distributions 
(Fig. 0]i-e) , one observes that the near-side peak topology remains essentially un- 
modified. We attribute this to a possible trigger bias of the near-side jet. However, 
it should be noted that the near-side yield does show a mild rise with centrality as 
discussed below. In contrast to the near-side peak, the away-side peaks show indi- 
cation for strong modifications (both in magnitude and in shape) for all centralities 
other than the most peripheral event selection. More precisely, in the 0-5%, 5-10% 
and 40-60% samples, the away-side peak evidences a plateau like shape which is 
decidedly non-gaussian and much broader in width than that for the 60-90% sam- 
ple. For the 10-20% and 20-40% samples (Fig. 0t-d), the away-side peak remains 
broad but also indicates an apparent local minimum at A(f> — 7r and a maximum 
at A</> = 27r/3. This latter pattern is similar to recent predictions of jet-induced 
"conical flow" [[23 1281 I29j . It should be pointed out however that these results 
do not preclude an alternative scenario which conjectures the combined influence of 
energy loss and the inclination angle of the jet with the flow field [122]. Nonetheless, 
both approaches require relatively strong coupling between jets and the high energy 
density medium. 

The solid (dashed) lines in Fig0] represent the conditional yield distributions 
that would result from subtracting out a v 2 product lowered (raised) by one system- 
atic error, respectively. The systematic error on v 2 is dominated by the uncertainty 
on the reaction plane dispersion. The dotted line indicates the jet-pair distribution 
resulting from a subtraction with v 2 product lowered by twice the systematic un- 
certainty From these curves, one can see that a decrease in v 2 by two intervals of 
the systematic error can recover the local minimum at Acf> = n. However, the result 
of a broadened away-side jet remains robust and the away-side jet-shapes remain 
non-Gaussian. Several studies are currently underway to firm up the mechanism/s 
responsible for the atypical away-side jet topologies observed in Au+Au collisions. 

To better quantify these jet properties, we split the jet-pair distributions into 
a near-side range, - A0 m ;„ and an away-side range A(f> m i„ - ir. The near- and 
away-side parts of the distribution are then further characterized by their RMS 
(taken around and tt) and their yield of associated pairs per trigger. These results 
are summarized in FigJSJas a function of centrality. For comparison, similar results 
are included for the 0-20% most central d+Au collisions (open circles) obtained at 
\J snn = 200GeV. The systematic uncertainty in v 2 has been propagated into the 
systematic errors for yields and RMS values. The systematic error on the yields 
also accounts for the systematic uncertainty on the single particle reconstruction 
efficiency. 

FigEb shows that both the near- and away-side widths for peripheral (60-90%) 
Au+Au collisions compare well with those obtained for d+Au collisions. This is 
consistent with the expectation of very little, if any, medium induced modifications 
to jet-topologies in peripheral Au+Au collisions. An inspection of the centrality de- 
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Fig. 5. (a) Associated yields for near- and away-side peaks in the jet-induced pair 
distribution; and (b) Widths (RMS) of the peaks. The open symbols denote results 
for 0-20% most central d+Au events from a recent PHENIX analysis. 



pendence of the near-side RMS shows no apparent change with centrality. On the 
other hand, the away-side width is significantly broadened for all but the most pe- 
ripheral event sample, possibly indicating strong modifications to the fragmentation 
process by the hot nuclear medium. Although the near-side widths are centrality 
independent, Fig|SJi points to a mild increase in the near-side conditional yield from 
peripheral to central Au+Au collisions, possibly indicating that even the near-side 
fragmentation process might still be influenced by the medium. The apparent dif- 
ferences in the evolution of near- and away-side jet characteristics could be signaling 
the contribution of several different mechanisms to jet-modification. 

5. Conclusions 

In summary, we have used a novel correlation function technique to decompose jet 
correlations from collective long range harmonic correlations (elliptic flow). We find, 
that the extracted jet-pair distributions show a strong centrality dependent change 
in shape and associated per-trigger yield, especially for the away-side jet. The jet- 
pair distributions obtained for peripheral Au+Au collisions are very similar to those 
obtained for d+Au collisions but the distributions for more central Au+Au collisions 
are markedly different and are qualitatively consistent with several recent theoretical 
predictions of possible modification to jet fragmentation by a strongly interacting 
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medium [ |2HI 03 HH1 HUH GUI- Further experimental and theoretical studies are 
clearly required to establish the detailed mechanism/s responsible for the observed 
jet modification/s and to pin down the properties of the high energy density strongly 
interacting matter produced in energetic Au+Au collisions at RHIC. Several such 
studies are currently being pursued with vigor. 
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